/***************************************************************************
 *   Copyright (C) 2006 by Harish Kukreja                                  *
 *   harish.kukreja@gmail.com                                              *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

#include <physiksbase/world.h>
#include <physiksbase/objects.h>

using namespace std;
using namespace PhysiksBase;

int main(int argc, char *argv[]) {
    World *world1 = new World(3, 1.0E-3);
//     world1->output = true;
    Vector *pos1 = new Vector(3, 0.0);
    Vector *pos2 = new Vector(3, 0.0);
    Vector *vel1 = new Vector(3, 0.0);
    Vector *vel2 = new Vector(3, 0.0);
    pos2->set(0, 6.3781E6);
    PhysicalObject *obj1 = new PhysicalObject(45.54, 0.0, *pos1, *vel1);
    PhysicalObject *obj2 = new PhysicalObject(5.9742E24, 0.0, *pos2, *vel2);
    world1->add_object(*obj1);
    world1->add_object(*obj2);
    world1->integrate(0.0, 5.0);

    delete pos1;
    delete pos2;
    delete world1;
    delete obj1;
    delete obj2;
    delete vel1;
    delete vel2;

    return 0;
}
